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We have investigated the doping dependence of the penetration depth vs. temperature in elec- 
tron doped Pr2- 3; Ce a ;Cu04_i using a model which assumes the uniform coexistence of (mean-field) 
antiferromagnetism and superconductivity. Despite the presence of a d x 2_ y 2 pairing gap in the 
underlying spectrum, we find nodeless behavior of the low-T penetration depth in underdoped case, 
in accord with experimental results. As doping increases, a linear-in-T behavior of the penetration 
depth, characteristic of d-wave pairing, emerges as the lower magnetic band crosses the Fermi level 
and creates a nodal Fermi surface pocket. 

PACS numbers: 74.20.Rp, 74.25. Dw, 74.20.Mn, 74.25. Nf 



An understanding of the symmetry of the order pa- 
rameter and its evolution with hole and electron doping 
is a key to unraveling the mechanism of high-T c super- 
conductivity in the cuprates. Many experimental and 
theoretical studies of these fascinating materials demon- 
strate the presence of antiferromagnetic (AFM) order in 
underdoping for both hole-[l[ and electron doping 0,0, H|. 
With hole-doping, the route followed by the AFM phase 
as it develops into the superconducting (SC) phase in- 
volves the intervention of nanoscale phase separations 
related to stripe or pseudogap physics The behav- 
ior with electron doping, on the other hand, seems to be 
simpler in that the doped phase appears to be a uniform 
AFM metal, possibly evolvinginto a phase with coexist- 
ing AFM and SC orders [1 SQ. 

For the hole doped cuprates it is generally be- 
lieved that d-wave pairing survives up to the edge 
of antiferromagnetism [J 0, @|, but the doping de- 
pendence of the pairing symmetry with electron dop- 
ing remains a matter of debate. This symmetry 
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sitive probes[l6|, in a variety of electron-doped cuprates, 
including Nd 2 _ a: Ce a ;Cu04_5 (NCCO)@, I i M 0, 
Laz-zCezCuCw (LCCO)^ [H Pr 2 _ a: Ce :c Cu04-5 
(PCCO)[H GS G3]. The results have been con- 
tradictory, with some early measurements 0, H, @, E3] 
finding evidence for s-wave pairing, while other exper- 
iments suggest a transition from d-wave in underdoping 



to either s-wave [111 1 1 31 ] or (d + is)-wave character [14j 



in the optimally and overdoped cases. Yet other 
experiments [Til. Il6| report only d-wave pairing, with the 
situation further complicated by the presence of non- 
monotonic SC-gap variations observed in NCCO[17| and 
Pri_ x LaCe x Cu0 4 -a (PLCCO)0. 

A recent study approximated the AFM background 
by treating the resulting partially-gapped Fermi surface 
(FS) in a two band model[3]. To understand the in- 
terplay between AFM and SC orders and the role of 



AFM order in modifying the pairing [19), |20|, in this arti- 
cle we directly evaluate the PD in a model with coexist- 
ing AFM and SC order. We assume a SC gap of d-wave 
pairing with a combination of first and third harmon- 
ics, which is necessary to incorporate nonmonotonic gap 
variations @, QjJ- We find that even in the presence of a 
d-wave pairing gap, the PD varies exponentially at low T 
for most dopings — a behavior characteristic of a node- 
less SC-gap, as antiferromagnetism suppresses the spec- 
tral weight from the nodal point. In the overdoped case 
(x = 0.152), the PD shows a linear-in-T behavior as the 
hole pocket forms in the nodal region. Our analysis in- 
dicates that with increasing electron doping the position 
of the maximum leading edge gap on the FS moves away 
from the antinodal direction and that the nonmonotonic 
nature of the gap becomes stronger. 

Our treatment of the in-plane PD is based on the 
Hamiltonian 



H — H.„, 



Hi, 



(1) 



where H pa i r describes the physics of coexisting AFM and 
SC orders. We take H va ir to be a one-band, tight-binding 



Hubbard Hamiltonian along the lines of Ref. [J in which 
the SC gap is of d-wave pairing with a combination of 
first and third harmonics. The tight binding parameters 
are assumed to be same as for NCCO|3|- The external 
perturbation is given by the electromagnetic interaction, 



ff « = -(-D^[E^ lB 



(2) 



where ct (c^ a ) is the electronic creation (destruction) 

operator with momentum fc, charge e and spin a, and 
c is speed of light. A is the Fourier component of the 
vector potential in momentum space, = d^/ihdk) is 
the band velocity for the noninteracting band 

The PD is obtained by evaluating the induced current 
parallel to the vector potential, which is proportional to 
the inverse square of the in-plane PD[2l|. We have gen- 
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eralized the pure BCS result to the mixed AFM-SC case and find[22 
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Here, a is the in-plane and d the out-of-plane lattice con- 
stant of PCCO and = l/k B T. The prime on the k 
summation means that the sum is restricted to wave vec- 
tors in the magnetic zone. The magnetic field is assumed 
to lie perpendicular to the CuC>2 plane. For a tetragonal 
lattice Xij is diagonal, with X aa = X^ = X within the 
CuC>2 plane. Interestingly, Eq. 3 displays a form simi- 
lar to that for a pure d-wave SCjH, excepting two 
modifications. Firstly, the FS has components v — ± as- 
sociated with the upper magnetic band (UMB) and the 
lower magnetic band (LMB), respectively: 
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where E Q% = 



sJ(^) 2 + (U Q S) 2 and e± = ± e £+ g)/2. 
is the SC gap and UqS the AFM gap in terms 
of the AFM repulsion Uq and the commensurate mag- 
netisation S at the nesting vector Q = Sec- 
ondly, the band masses and quasiparticle veloci- 
ties v"- have magnetic correlation corrections: 1/m- = 

k ' kij 

d2 (tt + uE^y^dkidkj), m = d($ + vE oi )/{hdk). 

We obtain the AFM and SC gaps self-consistently as a 
function of T for a series of doping levels over the range 
x = 0.115 — 0.152, using doping dependent interaction 
parameters^, before proceeding with the PD calculation 
from Eq. [3J The effective AFM interaction given by Uq 
is taken from our earlier work on NCCO[3j and decreases 
from a value of 3.33i at x=0.115 to 2>.lt at £=0.152; the 
resulting self-consistent magnetization S decreases lin- 
early from 0.2 to 0.13 over this doping range, in agree- 
ment with earlier results Q, despite the presence of the 
SC order. 

The two terms on the right hand side of Eq. [3] corre- 
spond to the conventional diamagnetic (first term) and 
paramagnetic 25| (second term) currents of electrons. In 
a London picture, A~ 2 (T) is proportional to the SC elec- 
tron density n s , and hence vanishes as T — > T c , while 
at T = all the electrons are superconducting. Here we 
find a similar result, but only the electrons in the AFM 
pockets condense. Similarly, as T —> 0, the linear-in-T 
PD found in overdoped samples reveals the presence of 
gap nodes, where normal quasiparticles persist to zero 
energy. 

Figure [T] compares the theoretical and experimental 
values of the inverse square of the PD in PCCO[26[ over 
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FIG. 1: (color online) Various colored lines give computed 
A -2 (T)/A -2 (0) as a function of T for different dopings x\ the 
corresponding experimental data for PCCO[26l] is shown by 
symbols of the same color (see legend). Inset: Computed 
(red line) and experimental (blue dots) values of A _2 (0) as a 
function of doping. Green arrow points to the kink associated 
with the opening of the nodal pocket in the theory curve. 



the doping range x = 0.115 — 0.152. The results are 
normalized to the computed T = value for the c-axis 
lattice constant d — 12.2 A 27] in order to highlight In- 
dependencies. The overall agreement is quite good, al- 
lowing us to adduce the doping dependence of the AFM 
and SC properties as discussed below. A discrepancy is 
found at the lowest and highest dopings, where the PD 
shows a tail extending beyond T c , possibly associated 
with sample inhomogeneities [28j| . 

Turning to the inset in Fig. 1, note first that the the- 
oretical values (red line) of A~ 2 (0) <x n s (0), do not in- 
volve any further fitting parameters beyond those used 
in fitting the T-dependence of A. Around optimal dop- 
ing, theory and experiment are seen to be in accord in- 
dicating that the theory correctly predicts the value of 
n s (0), although striking deviations are seen away from 
optimal doping. Insight into this behavior is obtained 
by observing that as T — * 0, all the electrons on the 
FS condense so that n s (0) is proportional to the area of 
the FS pockets. For this reason, the computed n s (0) 
decreases linearly with underdoping and undergoes a 
change in slope (marked by the green arrow) as the 
(7r/2,7r/2) pocket crosses the Fermi level in overdoping. 
In sharp contrast, the experimental points (blue dots) 



FIG. 2: (color online) (a) UMB contribution to the theoretical 
PD at three different dopings. Solid lines give the correspond- 
ing exponential fits of form [1 — Cexp(— (3D)] at low T. (b) 
LMB contribution at x — 0.152 (circles), and the related lin- 
ear fit, (1 — fflT), at low T (blue line). 



present a peak around optimal doping and a loss of car- 
riers away from optimal doping, indicative of 'bad metal' 
physics where the SC transition is dominated by thermal 
phase fluctuations 2^, 3(| • These results suggest that on 
the underdoped side AFM fluctuations are more dele- 
terious than expected from the mean field BCS model 
underlying our computations. The reappearance of bad 
metal behavior on overdoping is puzzling and its origin is 
unclear— it may be related to increasing doping-induced 
disorder. 



We discuss the doping and T-dependencies of the pre- 
ceding theoretical PD results with reference to Figures 
[2]4. In Fig. Ufa), we emphasize that the contribution of 
the UMB is nodeless since the PD is dominated by ener- 
gies near the Fermi level and the UMB pocket is far from 
the nodal region. An exponential form, 1 — C exp (—/3-D), 
is seen to produce an excellent fit in the low-T region in 
Fig. da). The values of the SC-gap D in Fig. E(a) 
so obtained for the UMB (blue dots) are quite close to 
the SC-gap (red open squares) at the electron pocket tip, 
marked by yellow diamonds on the FS plots of Figs.BJb)- 
(d). The value of C is «4.4, essentially independent of 
doping. Only at the highest doping x = 0.152 do we find 
a significant linear-in-T (d-wave) contribution to the PD 
as shown in Fig. HJb), which coincides with the appear- 
ance of the hole pocket near (tt/2,w/2) at high dopings as 
the LMB crosses the Fermi level (see Fig. 4(d)). A linear 
equation of the form, 1 — aT, fits the LMB contribution 
very well up to T — 1.5K as shown by blue line in Fig. 
HJb) [3l[ . We do not find a second regime of linear-in-T 
PD in the strongly underdoped regime [Tij. 

Figure [3] examines the doping dependence of the SC 
gap parameters. The dome-like shape as a function of 
doping of the gap D in Fig. 02a) is reflected in the be- 
haviors of the first and third harmonics of the pairing 
gap in Fig. [3](b) as well. Fig. 3(c) delineates the dop- 
ing dependence of the first and third harmonics of the 
pairing interaction, which display a maximum near x ks 
0.13 where T c is optimal. The doping dependence of the 
SC gap parameters Ai and A3 is compared with vari- 
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FIG. 3: (color online) (a) SC-gap D as a function of doping 
obtained from exponential fits similar to those in Fig. 2(a) 
at low T (blue dots), and the values of the SC gap (red 
open squares) at the tip of the electron pockets, shown in 
Fig. |4jb)-(d) (yellow diamonds), (b) Self-consistent first har- 
monic 2Ai (blue dots) and third harmonic — 4A3 (red open 
circles) of the SC gap are compared with several experimental 
results for NCCO and PCCO: Raman scattering in B29 chan- 
nel from NCCO [PCCO] [321L green filled [green open] squares; 
Tunnelling data on NCCO[13l. Il4jl. magenta diamonds; and 
two band model computations 18], orange triangles. Green 
dashed line shows the scaled Raman B29 channel gap for 
NCCO. (c) SC interaction potentials as a function of doping: 
Absolute value of the d-wave first harmonic — Vi (blue dots) 
and the third harmonics 2.3V3 (red circles), (d) 2A/fcsT c for 
Ai and — I.5A3 for the first and third harmonic SC gaps are 
compared to the experimental results. Various symbols have 
the same meanings as in (b). 



ous experimental results in Figs. 3(b) and (d). Some 
disagreement with Raman experiments on NCCO (green 
filled squares) and PCCO (green open squares)^ is due 
to sample variations, reflected in T c variations, while the 
ratio 2Ai/ksT c is essentially constant and agrees well 
with experiment. If we scale the experimental gap to fit 
the calculated maximum at optimal doping, we can re- 
produce the dome-like behavior of the SC gap as shown 
by the green dashed line in Fig [3{b) for NCCO. The 
tunnelling data (magenta diamonds) [lj, [3| do not show 
a maximum, because tunnelling is sensitive to the total 
gap obtained by combining AFM and SC gaps, and this 
combined gap in our computations does not have a maxi- 
mum near optimal doping. Similarly, the larger gap seen 
by Raman [32| in under- and optimally doped samples 
can be understood since the Bi g channel measures the 
total spectral gap near the (tt, 0) point (33|, and hence is 
strongly coupled to the AFM order. 

Fig. Eld) shows that the ratio 2Ai//cbT c possesses a 
nearly constant value of 4.1, close to the BCS value for 
a d-wave gap. In contrast, for fixed ratio of Vi/Va[34j|. 
the third harmonic ratio — 2A3/fcsT c increases linearly 
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FIG. 4: (color online) a) Variation in the leading edge gap 
(LEG), A s hift, over the FS in terms of the FS angle <fr, where 
ij> increases from zero along the antinodal direction to 45° 
along the nodal direction. Solid dots mark positions of the 
LEG maximum, </> ma x, which is seen in the inset to increase 
nearly linearly with doping, (b)-(d) FS is calculated for three 
different dopings x. Red color denotes maximum spectral 
weight and blue color zero intensity. Red dots give the posi- 
tion of the maximum LEG gap on the FS considered in (a). 
Yellow diamonds mark the tips of the UMB electron pockets, 
at which the SC gaps in Fig. 2(c) were calculated. 



with doping. This is the reason that the position of the 
maximum of the leading edge gap (LEG) A s hift on the 
FS, given by the FS angle </> inax hi Fig. 4(a), moves away 
from the antinodal point with doping; interestingly, the 
hot spots also move away with doping from the antin- 
odal direction, but their shift is much smaller. We find 
that the ratio of the maximum value of A s hift to its value 
along the antinodal direction increases with doping, indi- 
cating that the non-monotonic nature of d-wave pairing 
symmetry becomes more pronounced as one goes from 
under- to overdoping in the electron doped cuprates. 

The evolution of the (tt/2, 7r/2)-centered nodal holc- 
pocket is seen in Figs. Hfb)-(d). The absence of nodal 
pockets in the underdoped regime (see (b)) is responsi- 
ble for the nodeless behavior of the SC gap. At optimal 
doping x=0.137, the hole pocket is still ~ 25 meV be- 
low Ep, but can be seen in (c) due to the finite energy 
resolution. The nodal pocket is fully formed in the over- 
doped case of (d) which is related to the striking c?-wave 
behavior of PD in Fig. 2(b) as well as the kink in n s (Q) 
in Fig. 1 inset. 



In conclusion, we have shown that the linear-in-T vari- 
ation of A~ 2 in electron doped cuprates is related to the 
appearance of the (7r/2, 7r/2)-nodal hole pocket on the 
FS, which occurs in the overdoped regime. In under- 
doping, where the FS only consists of the (n, 0)-centered 
electron pockets, A~ 2 varies in a nodeless manner, even 
though the pairing interaction is of <i-wave symmetry, 
because the electron pocket lies far from the nodal re- 
gion. Our analysis indicates that the SC electron density 
(n s (0)) is suppressed in a non-BCS fashion as one goes 
away from optimal doping to either under- or overdoping. 
Interestingly, we find that the SC interaction (Vi and V3) 
also peaks at optimal doping. 
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